Event-by-event charged-neutral fluctuations in Pb-i-Pb collisions at 158 A GeV 

(WA98 Collaboration) 



M.M. Aggai-wal", Z. Ahammed'', A.L.S. Angelis"^'', V. Antonenko'', V. Arefiev'^, V. Astakhov% V. Avdeitchikov'', T.C. Awes^ 
P.V.K.S. Babas, S.K. BadyaF, S. Bathe\ B. Batiounia^ T. Bemier', K.B. Bhallai, VS. Bhatia\ C. Blume\ D. Bucher^ 
H. Busching'', L. Carlen'', S. Chattopadhyay'', M.P. Decowski', H. Delagrange', P. Donni'', M.R. Dutta Majumdai"'', 
K. El Chenawi'', A.K. Dubey'', K. Enosawa", S. Fokin'', V. Frolov=, M.S. Ganti'', S. Garpman'' ', O. Gavl^shchuk^ RJ.M. Geurts°, 

T.K. Ghosh'', R. Glasow'', R. Gupta^, B. Guskov^ H. A.Gustafsson'^'i, H. H.Gutbrodi, I. Hrivnacova', M. Ippolitov'', 
H. Kalechofsky"^, R. Kamermans°, K. Karadjev'*, K. Karpio\ B. W. Kolb'', I. Kosarev*^, I. Koutcheryaev'^, A. Kugler', P. KuHnich', 
M. Kurata", A. Lebedev'', H. Liu\ H. LohnerP, L. Luquin', D.P Mahapatra", V. Manko'', M. Martin^ G. Martinez', A. Maximov", 
Y. Miake", G.C. Mishra™, B. Mohanty'', M.-J. Mora', D. Morrison', T. Moukhanova'', D. S. Mukhopadhyay^ H. Naef, 
. . B. K. NandP, S. K. Nayaks, T. K. Nayak'', A. Nianine'', V. Nikitine=, S. Nikolaev'', P Nilsson'', S. Nishimura", P Nomokonov^ 
^""^ J. Nystrand*^, A. Oskarsson'', I. Otterlund'^, S. PavUouk'', T. Peitzmann", D. Peressounko'', V. Petracek', V. Petracek"", 

^ '. W. Pinanaud', R Plasil^, M.L. Purschke'', J. Rak^ R. RaniwaW, S. RaniwaW, N.K. Rao^, E Retiere', K. Reygers'', G. Roland', 
^ L. Rosselef , I. Roufanov^ C. Roy', J.M. Rubio^ S.S. SambyaF, R. Santo'', S. Sato", H. Schlagheck'', H.-R. Schmidt^, Y. Schutz', 
G. Shabratova^ T.H. Shah^, I. Sibhiak'', T. Siemiarczuk^ D. Silvemyr'', B.C. Sinha'', N. Slavine^ K. Soderstrom'', G. Sood% 
S.P. S0rensen', P. Stankus^, G. Stefanek", P. Steinberg', E. Stenlund*", M. Sumbera^ T. Svensson'', A. Tsvetkov'', L. Tykai'ski", 
E.C.v.d. Pijll°, N.v. Eijndhoven°, G.J. v. Nieuwenhuizen', A. Vinogradov'', YP Viyogi*", A. Vodopianov'^, S. Voros'^, B. Wyslouch', 

G.R. Youngf 



X 



> 
00 

(N 

rn 
O 



"University ofPanjab, Chandigarh 160014, India 
^Variable Energy Cyclotron Centre, Calcutta 700 064, India 
University of Geneva, CH-I2I1 Geneva 4, Switzerland 
''RRC Kurchatov Institute, RU-123182 Moscow, Russia 
" Joint Institute for Nuclear Research, RU-141980 Dubna, Russia 
f Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6372, USA 
^University of Jammu, Jammu 180001, India 
^ University ofMUnster, D-48I49 Munster, Germany 
'SUBATECH, Ecole des Mines, Nantes, France 
'University of Rajasthan, Jaipur 302004, Rajasthan, India 
''Lund University, SE-221 00 Lund, Sweden 
'MIT Cambridge, MA 02139, USA 
"'Institute of Physics, 751-005 Bhubaneswar, India 
" University ofTsukuba, Ibaraki 305, Japan 
"Universiteit Utrecht/NIKHEF, NL-3508 TA Utrecht, The Netherlands 
''KVI, University of Groningen, NL-9747 AA Groningen, The Netherlands 
1 Gesellscliaft fiir Schwerionenforschung (GSI), D-64220 Darmstadt, Germany 
''Nuclear Physics Institute, CZ-250 68 Rez, Czech Rep. 
' Soltan Institute for Nuclear Studies, PL-00-681 Warsaw, Poland 
' University of Tennessee, Knoxville, Tennessee 37966, USA 



X 



Abstract 

.' Charged particles and photons have been measured in central Pb+Pb collisions at 158 A GeV in a common (// - 0)-phase space 
region in the WA98 experiment at the CERN SPS. The measured distributions have been analyzed to quantify the frequency with 
which phase space regions of varying sizes have either small or large neutral pion fraction. The measured results are compared with 
VENUS model simulated events and with mixed events. Events with both large and small charged-neutral fluctuations are observed 
to occur more frequently than expected statistically, as deduced from mixed events, or as predicted by model simulations, with the 
difference becoming more prominent with decreasing size of the Arj - region. 
PACS: 25.75. -q,25.75.Gz 

Keywords: The WA98 experiment. Photon Multiplicity Detector, Disoriented Chiral Condensate, Silicon Pad Multiplicity 
Detector,Sliding Window Method, Charged-Neutral Fluctuations. 
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Enhanced fluctuation in the production of neutral versus 
charged pions has been one of the predicted signals for chiral 
symmetry restoration in heavy ion collisions at ultra-relativistic 
energies UjJ. It has been proposed that the extreme energy 
density of the matter produced in the region of spatial over- 
lap between two heavy ions colliding at relativistic energies 
may provide the physical conditions necessary for the forma- 
tion of a chiral condensate that may be aligned in a direc- 
tion different from the true vacuum. Domains of such Dis- 
oriented Chiral Condensates (DCCs) are expected to emit pi- 
ons coherently from the collision volume, which may result in 
large fluctuations in the neutral pion fraction, /, defined as / 
= Njr"/iN„i> + N„±), where A^;^" and A^;^* are the multiplicities of 
the neutral and charged pions, respectively. The neutral pion 
fraction / for DCC domains is predicted to follow a probabil- 
ity distribution of the form P(/) = 1/2 ^Jf |2], which is very 
different from the binomial (or negative binomial) distribution 
observed for generic pion production in hadron collisions. 

Events with large charged-neutral fluctuations, the so called 
Centauro (Nd, » Ny) and Anti-Centauro (Ny » Nch) events, 
were first observed in the JACEE cosmic ray experiment [jsl. 
The search for such unusual events and DCC formation was 
later carried out by the DO fl, CDF fl, and Minimax fl ex- 
periments in p-Hp collisions at the Tevatron and by the NA49 
m and WA98 M, 

ll iL [121] experiments in nuclear collisions 
at die CERN SPS. Upper hmits to flie production of DCC do- 
mains have been set within various model assumptions. WA98 
has conducted an exhaustive search for charged-neutral fluctu- 
ations in nuclear colUsions. In the first WA98 analysis ifiol . 
correlated neutral vs. charged particle fluctuations were inves- 
tigated globally over a large tj - (p acceptance of more than one 
unit of rapidity near mid-rapidity, with full azimuthal coverage. 
Subsequent investigations based on statistical correlations yj | 
or multi-resolution Discrete Wavelet Transform (DWT) tech- 
niques II12II . searched for DCC-like fluctuations in localized 
regions of {rf - <p) phase space. 

In this article we present an improved search for localized 
charged-neutral fluctuations in central Pb + Pb collisions at the 
SPS using a simple SHding Window Method (SWM) 1 13]. This 
method utilizes the full azimuthal resolution of the WA98 detec- 
tors to identify all regions with unusually large or small values 
of the neutral particle fraction. The sensitivity of the SWM has 
been studied using simulated data based on simple model as- 
sumptions of charged-neutral fluctuations similar to those used 
in the study of the DWT method illllli] . The SWM is found to 
have a sensitivity which is limited only by the available statis- 
tics and can be several orders of magnitude better than the DWT 
method for a given data sample. Preliminary results have been 
presented in Refs. 1 15, 16.1 . 

The analysis reported here used the measurement of pho- 
tons and charged particles produced in Pb-i-Pb collisions at 158 
A GeV recorded by the WA98 experiment, during a period of 
magnetic field-off' operation of WA98 in the 1996 run period of 
the CERN SPS. Photons were measured with the Photon Mul- 
tiplicity Detector (PMD) [ 17] and charged particles were mea- 
sured with the Sihcon Pad Multiphcity Detector (SPMD) id]. 
The analysis is restricted to their common region of rj - (p ac- 



ceptance overlap. The PMD was located at 21.5 meters down- 
stream of the target and consisted of plastic scintillator pads 
of varying sizes arranged inside 22 box modules placed behind 
3X0 thick lead converter plates It covered the pseudo- 

rapidity region 2.9 < 77 < 4.2. For the photon identification 
criteria applied the average photon counting efficiency and the 
purity of the photon sample were found to be 68% and 65%, 
respectively, for the central event sample used lll7ll . The az- 
imuthal resolution of the PMD is much less than 1° and the limit 
for detection of low pj particles is 30 Me V/c. The SPMD was 
located at 32.8 cm from the target and consisted of 22 radial and 
46 azimuthal segments in each of four quadrants. It covered the 
pseudorapidity region of 2.35 < 7; < 3.75 ]JJ. The charged par- 
ticle detection efficiency of the SPMD was 99%. The azimuthal 
resolution of the SPMD is 2° and the limit for detection of low 
Pt particles is nearly zero. The total transverse energy Ej mea- 
sured in the pseudorapidity region 3.5 < 77 < 5.5 by the Midra- 
pidity Calorimeter (MIRAC) jfisl], located 24 m downstream 
of the target, was used to characterize the event centrality. All 
procedures for event selection and removal of backgrounds de- 
scribed in detail in previous publications I^ljlll2lll5lll7|] were 
followed in the present analysis. A total of 185 K events, be- 
longing to the 15% most central collisions of the WA98 mini- 
mum bias cross section, have been analyzed. 

In the present analysis, the multiplicity of photons measured 
in the PMD and the multiplicity of charged particles mea- 
sured in the SPMD, in the region of common 77 acceptance 
(2.9 < 7/ < 3.75), are employed as experimental observables 
to approximate the neutral pion fraction by the neutral particle 



fraction / 



Nyl2 



The neutral particle fraction / is calcu- 



Ny/2+N,i, ■ 

lated within windows of varying size A<p in azimuth. 

The search for non-statistical fluctuations requires compari- 
son with a model-independent baseline, free of dynamical cor- 
relations, that includes detector effects and fluctuations of sta- 
tistical nature only. For the analysis of experimental data the 
technique of event-mixing can provide such a reference data 
sample in which dynamical correlations among particles are 
completely destroyed, leaving only correlations that might re- 
sult from efficiency variations. The mixed events are generated 
by the standard procedure of using real events to construct new 
mixed events with particles selected randomly, one particle per 
real event. In this analysis, the mixed events are created with 
the same global multiplicity distribution as the experimental 
data. Furthermore, for each real event, a mixed event is cre- 
ated with the same global charged particle and photon hit mul- 
tiplicity, thereby maintaining the same global charged-neutral 
multiplicity correlation. Two kinds of mixed events are used in 
this study 112]. In the first case, termed Ml mixed events, all of 
the charged particle and neutral hits in each mixed event are se- 
lected from different real events. Thus the Ml mixed events 
constitute fully mixed events with only the global charged- 
neutral multiplicity correlation maintained. The second class of 
artificial events, termed M2 mixed events fl^, are created by 
mixing PMD and SPMD subevents from different events, but 
without mixing hits within the individual detectors. In the M2 
mixed events the charged and neutral fluctuations are separately 
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Figure 1: Top: The / distributions for data and mixed events for tlie 5% most 
central Pb+Pb events for randomly selected rj- </> regions with A<l> =60° . A fit to 
the distribution for mixed events with a Gaussian function, with fit parameters 
mean = 0.33923±0.00004 and cr =0.04041±0.00004 , is also shown. Bottom: 
Ratio of the / distributions for the data and the Gaussian fit to the mixed events 
(open circles) and the ratio of the mixed events with its Gaussian fit (solid cir- 
cles). 

maintained, but the correlation between charged and neutral hits 
is removed, except that the global charged-neutral multiplicity 
correlation is maintained. In what follows the mixed events are 
the Ml type unless explicitly stated. 

Two approaches have been used to characterize the statistical 
fluctuations and trivial charged-neutral fluctuations arising in 
generic particle production. Statistical fluctuations due to finite 
particle effects and small data samples are studied by a compar- 
ison with mixed events generated from the data, as mentioned 
above. In addition, the magnitude of trivial charged-neutral 
fluctuations arising in generic particle production have also 
been studied with simulated events produced by the VENUS 
event generator [ 19] passed through the GEANT detector sim- 
ulation fcoil of the WA98 experiment setup, here referred to as 
VG events, and mixed events generated from these VG events. 
For this study, 27k VG simulated Pb-t-Pb collision events have 
been generated with impact parameter range selection corre- 
sponding to the 15% most central collisions of the data sample. 

Fig.[T](top) shows the / distribution for data for the 5% most 
central Pb-i-Pb events within different regions of A(p -60° in 
an event. For comparison, the distribution obtained from Ml 
mixed events is also shown together with a Gaussian function 
fitted to the mixed event distribution. It is observed that the / 
distributions for the data and mixed events are asymmetric and 
extend beyond the Gaussian fit to larger /-values. In Fig. [1] 
(bottom) the ratio of the / distribution for the data to the Gaus- 
sian fit to the mixed event distribution is shown along with the 
ratio of the mixed event distribution to its Gaussian fit. The 
high value of the ratio at large / signifies vanishing small val- 
ues of the fitted distribution. The / distribution for the data is 



Figure 2: The fi„a_x (left) and /„„„ (right) distributions for A<f) = 60° window 
size for data (top) and VG (bottom) for the 5% most central events together 
with the distributions for their respective mixed events. 

observed to be broader than that of the mixed events with small 
and large /-values occurring more frequently. Because of the 
predominance of generic pion production around / ^ 1/3, non- 
statistical fluctuations can best be studied in the tail regions of 
the / distribution. 

For this reason we focus our attention on the maximum and 
the minimum /-values (i.e., f,„ax and /,„>,) amongst all Arj-Acp 
regions in each event in order to search for photon-excess and 
charge-excess type fluctuations. This is done using the SWM, 
where the azimuthal plane is scanned by sliding the A(p window 
of chosen size in steps of 6(p - 2° (limited by the azimuthal 
resolution of the SPMD) and computing / for each A0 region 
to extract the maximum and the minimum values of / in each 
event, represented by f„,ax and /,„„, respectively. The window 
size is kept fixed in rj and varied in azimuth from A0 - 20° 
to 150°. The f,„ax and /,„>, distributions for the experimental 
data and for VG events for the 5% most central Pb-nPb events 
are shown in Fig. |2]for a window size of A(p - 60°. The dis- 
tributions for mixed events generated from data and VG events 
are also shown. The measured /„„( and /,„■„ distributions are 
seen to extend further to the right and left, respectively, than the 
mixed events distributions, with the difference observed in the 
data seen to be greater than that for the VG events. 

In order to study the fluctuations in greater detail the mean 
(fj.) and RMS deviation (cr) are calculated for the mixed event 
distributions for both the experimental data and the VG events. 
The events in the tails of the f„ax and /,„„ distributions of the 
experimental data and the VG events are selected by applying 
a cut at yU ± ncr, with positive sign being appUed for the f^ax 
distribution and negative sign for the /„„„ distribution. Events 
having regions with f„,ax or /,„■„ values beyond the cut are la- 
beled as 'exotic' events. Results for n values of 3 and 4 are 
presented. 
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Figure 3: Scatter plot of Ny versus N^/, for a randomly selected 60° region for 
the 5% most central Pb+Pb events (solid points). The Ny versus Nrf, coiTelation 
for the photon and charge excess regions found by the SWM with the 3a" cut 
for the same event sample is also shown (open points). 



Figure 4: The percentage of events having A;;-A0 regions with fi„ax exceeding 
a 3iT cut (left) and with less than a 3it cut (right) versus the window size 
for the 5% most central events for data (top) and for VG (bottom) together with 
the results for their respective mixed events. Errors bars denote both statistical 
and systematic contributions. 



Fig. |3]displays a scatter plot of Ny versus Nc/, for a randomly 
selected = 60° region in each event. A weak correlation 
of Nc/, with Ny is observed due to residual impact parameter 
correlations, but with a wide spread. It is observed that a few 
points lie beyond the main cluster of events corresponding to 
events with large fluctuations in the multiplicities of photons 
and charged particles. When the SWM is applied to the same 
events to select those events where the / values exceed + 
3cr) for fmax or fall below (yu - 3cr) for /„„■„ many more events 
with regions having large fluctuations in the number of photons 
and charged particles are found, as shown by the open points in 
Fig.E] 

Fig. m shows the percentage of events beyond the /i + 3cr 
cut as a function of the window size A0. Results are shown 
for data and mixed events in the top panel and for VG events 
and VG mixed events in the bottom panel. The results for the 
mixed events from data and VG are seen to be very similar. If 
the f,„ax ifmin) distributions were well described by a Gaussian 
function, the percentage of mixed events exceeding the 3cr cut 
is expected to be 0.135%, independent of A(p. Instead, it is ob- 
served that due to the asymmetry of the /„,ni (/„„„) distribution 
(see Fig. |2]l there are excess number of events in the region 
of larger (smaller) /-values, with an excess that grows with 
decreasing A0 bin size, even for mixed events. It is observed 
that the percentage of exotic events in the data exceeds that in 
mixed events for all A(^ bin sizes for both the photon-excess 
and charge-excess selections. While the mixed events for data 
and VG events are very similar, the difi'erence between data and 
mixed events exceeds the difiference between VG and its mixed 
events. It is also seen that the difference between VG and its 
mixed events is much less for the charge-excess case than for 
the photon-excess case. The PMD-SPMD detector combination 



is better suited to study the photon-excess fluctuations because 
the purity of the photon sample in the PMD improves for large 
/ and can reach values up to 85% |,21]. For the charge-excess 
case the charged particle contamination plays a significant role 
in reducing the purity of the photon sample, which can be as 
low as 45%, and also introduces unwanted correlations due to 
charged hits in the SPMD that are simultaneously mis-identified 
as photon hits in the PMD. 

The various sources of systematic errors associated with A^y 
and Nch distributions have been investigated and described in 
detail previously [17]. The systematic error in the determina- 
tion of photon multiplicity is (-7.1%, -1-3.4%) and in the charged 
multiplicity it is +4%. To investigate the systematic errors on 
the number of observed patches new set of event samples was 
generated by randomly (a) removing 7.1% of photon hits and 
adding 4% of charged particle hits and (b) adding 3.4% photon 
hits and removing 4% charged particle hits. For these samples 
mixed events samples were also generated. These new event 
samples were analyzed in the same way as the real data and 
mixed events samples mentioned earlier to obtain the system- 
atic errors. Errors shown in figures are quadratic sum of statis- 
tical and systematic errors. 

The percentages of events having exotic photon and charge 
excess regions have been further investigated for additional 
types of event samples. In order to determine if the larger fluc- 
tuations observed in the data might be due to elliptic flow, the 
mixed and VG events have been modified to introduce ellip- 
tic flow following the method of Poskanzer and Voloshin 11221] . 
The pr dependent values of V2 are taken from Ref. fl^. Since 
the Pt- of the particles is not measured in the PMD or SPMD, a 
V2 value has been assigned to a given particle following the pj- 
probability distribution obtained from the VENUS event gener- 
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Figure 5: Percentage of events having regions witli fi„a_x (left) exceeding 3ct 
(top) or 4(T (bottom) cuts and with f„,i„ (right) less than 3cr (top) or Act (bottom) 
cuts versus Ad) for different sample types as desciibed in the text. 



Table 1 : Percentages of exotic events for the 4a- cut on the fi„„x (photon excess) 
or fiiiii, (charge excess) for Data, Ml mixed, Venus + Geant (VG) and M2 mixed 
events for different centralities for the A0 = 60° window size. The statistical 
and systematic errors are indicated. 



Cent. 


Sample 


fmax 


fmin 


0-5% 


DATA 


0.316±0.020 


0.105+0.012 




SYS ER 


+0.005-0.015 


+0.019-0.000 




Ml 


0.081+0.010 


0.014+0.004 




VG 


0.100+0.027 


0.028+0.014 




M2 


0.471+0.027 


0.083+0.011 


5-10% 


DATA 


0.260+0.019 


0.055+0.009 




SYS ER 


+0.016-0.006 


+0.011-0.003 




Ml 


0.048+0.008 


0.013+0.004 




VG 


0.099+0.033 


0.000+0.000 




M2 


0.471+0.028 


0.113+0.014 


10-15% 


DATA 


0.274+0.028 


0.045+0.011 




SYS ER 


+0.023-0.020 


+0.00-0.014 




Ml 


0.051+0.0140 


0.014+0.007 




VG 


0.074+0.033 


0.015+0.015 




M2 


0.357+0.032 


0.083+0.016 



ator. Fig. |5] (top) shows the percentage of events in which f^ax 
(left) exceeds the 3cr cut or in which /,„,„ (right) is less than 
the 3cr cut versus for Ml mixed events with flow (MI+V2) 
and VENUS+GEANT wifli flow (VG+V2). It is seen fliat flie 
observed enhancement in the percentage of exotic events in the 
data can not be explained as due to elliptic flow. 

The analysis was also performed on M2 mixed events con- 
structed as described above 111 ll 11211 by mixing the unaltered 
PMD hits of one event with the unaltered SPMD hits of a dif- 
ferent event, keeping the same total number of photons and 
charged particles as in the original event, and hence keeping the 
overall correlation between them. The results for the M2 event 
analysis shown in Fig. |5] are seen to lie above the data signify- 
ing that the charged-neutral fluctuations in the data tend to be 
correlated and give rise less frequently to events with photon- 
excess or charge-excess regions than the randomly correlated 
charged and neutral hits of M2 mixed events. This behaviour 
was previously reported in Ref. and is contrary to the naive 
expectations for a DCC. The result might be attributed to resid- 
ual impact parameter correlations in the charged and neutral 
multiplicity in the data. 

The curves shown in Fig.|5]were obtained from VENUS sim- 
ulations without GEANT detector response for the 5% most 
central collisions corresponding to the impact parameter range 
0-3.5 fm. The pure VENUS events are seen to follow closely 
the Ml mixed event results indicating essentially no correla- 
tions in the particle production in the VENUS event genera- 
tor The increase in the percentage of events with neutral or 
charge excess after application of the detector response (VG) 
is attributed to correlation due to detector effects where some 
charged particles may be mis-labeled as photons in the PMD. 

The bottom panels of Fig. |5] show the same analysis as in 
the upper panels but with the cuts on the f^^ix (left) and /j„,„ 



(right) distributions increased to 4cr. Although the percentage 
of events in the photon-excess or charge-excess regions are de- 
creased with the more stringent cuts, as expected, the general 
trends and conclusions are the same when comparing the re- 
sults for the 3cr cuts (Fig. |5](top)) with those for the 4cr cuts 
(Fig.|5] (bottom)). 

The analysis has also been performed for events in the 5-10% 
and 10-15% most central event selections. Table [1] compares 
these results for the Act cut with the 0-5% results for differ- 
ent event sample types as discussed above. This comparison 
shows little centrality dependence for the results over this lim- 
ited range of centralities, there being a slight tendency of in- 
creased excess in the data with increasing centrality. 

In addition to the event selection and data clean-up methods 
that have been employed, a set of additional checks have been 
performed with the 5% most central collisions for the = 60° 
window size to further examine the event structure and to search 
for possible detector artifacts. 

The first check concerns the azimuthal distribution of the ex- 
otic regions. Fig. |6] (top) shows the azimuthal distribution of 
the center of the azimuthal window observed to contain an ex- 
otic photon or charge excess region. The regions are observed 
to be distributed throughout the full azimuth, indicating that the 
exotic regions are observed throughout the entire acceptance 
regions of the detectors. 

The second check concerns the investigation of possible 
time-specific detector readout malfunctions. For each event 
found to have an exotic region, the neutral particle fraction (/) 
was calculated in the immediately preceding and immediately 
succeeding events in the same (Ari,A(f)) region. The / distribu- 
tions of these events associated to corresponding exotic events 
with /„,fl v values exceeding the 3cr cut or /„„■„ values less than 
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Figure 6: Results of various checks: (top) azimuthal distributions of exotic 
plioton and charge excess regions; (center) / distributions for preceding and 
succeeding events for the same region size at the same azimuth as indicated; 
and (bottom) charged particle multiplicity distribution in the non-overlapping 
region of the SPMD at the same azimuthal position as the exotic region and in 
a randomly selected region of the SPMD of similar size in the same event. 



the 3cr cut are shown in Fig. |6] (center). It is observed that the 
associated / distributions have their peak around 1/3 and ex- 
hibit behaviour similar to that of generic events, which in turn 
suggests normal detector operation. 

The third check investigated the operation of the SPMD by 
examining if the observed variation in the number of charged 
particles in an exotic region extended to the pseudorapidity re- 
gion 2.35 < T] < 2.9 of the SPMD that does not overlap with 
the PMD. Fig. |6] (bottom) compares the charged particle multi- 
plicity distributions in the non-overlapping part of the SPMD in 
the same azimuthal zone as the exotic region, and in the same 
event, with the distributions in a randomly selected region of 
similar size, for photon excess as well as for charge excess re- 
gions. The good agreement of the charged particle multiplicity 
distributions for regions near the exotic location and random re- 
gions suggests normal operation of the SPMD when the exotic 
regions are found. 

In summary, photons and charged particles have been mea- 
sured in a common 77 - phase space region in the WA98 exper- 
iment for 158 A GeV Pb+Pb collisions. The results were ana- 
lyzed for the 15% most central collisions to search for Arj-Aif) 
regions of photon or charged multiplicity excess using the slid- 
ing window method. The neutral particle fraction (/ ) has been 
calculated in Ari-Acp regions within 2.9< t] < 3.75 for different 
Aip windows. It is found that regions with photon or charged 
multiplicity excess occur more frequently in the data than in 
fully mixed events or in VENUS -hGEANT simulation events. 
The percentage of events with photon-excess or charge-excess 
increases with decreasing A<p bin size, but varies little if at all 
as a function of centrality within the 15% most central colli- 
sion event sample. The exotic regions having large charged- 



neutral fluctuations are distributed throughout the azimuthal ac- 
ceptance. Several investigations have been performed to rule 
out possible detector artifacts. The comparison with the anal- 
ysis of mixed events provides a model independent demonstra- 
tion of non-statistical charged-neutral fluctuations. 
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